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SUMMARY 


Inspection of the directional scintillation counter data from Explorer IV re- 
vealed fluctuations in the counting rate that appeared to be due to the body motions of 
the satellite and the anisotropy of the radiation flux. It was felt that the body motions 
and the angular dependence of the radiation flux could be determined from the data 
available. 

This work is a description of the methods used and results obtained from this 
determination. It was found that the counting rate is generally a maximum when the 
detector axis is normal to the magnetic field which is expected from the theory of 
particle trapping. The counting rate falls off as the detector axis makes smaller 
angles with the magnetic field lines, reaching a minimum value when the detector 
comes within about 40 degrees of the magnetic field lines. Knowledge of this angular 
dependence of the counting rate will allow calculation of the radiation distribution 
along a particular magnetic line of force. 

In the analysis of the rigid body motion, two rather surprising facts were ob- 
served. First, the components of the total angular momentum vector of the satellite 
do not remain constant in a space orientation-fixed reference system, although the 
magnitude remains essentially constant. Second, the small residual roll of the satel- 
lite about its longitudinal axis does not decay monotonically. In fact, at one point, the 
roll rate was observed to increase by almost an order of magnitude in one day. These 
two observations suggest that external forces are present and exert body torques on the 
satellite. It is possible that the origin of the forces may be due to an interaction of a 
magnetic moment in the satellite with the geomagnetic field, although this is not yet 
confirmed. 



SECTION I.. INTRODUCTION 


One of the most significant contributions from the satellite and space program 
during the International Geophysical Year was the discovery of the two radiation belts 
encircling the earth. While much has been learned about the configuration, compo- 
sition, and radiation intensities of these belts, comparatively little is known about the 
angular distribution of particle flux about magnetic lines of force. It is known that the 
radiation exhibits a strong anisotropy, but a detailed knowledge of intensity as a func- 
tion of angle relative to geomagnetic field lines will aid considerably in interpreting 
some of the observed results and in mapping the radiation belt. 

Satellite 1958 Epsilon, launched on July 26, 1958, contained two highly direc- 
tional scintillation detectors as part of its radiation experiment. The rigid body 
motion of this satellite was such that the detector axes traced out complex motions on 
a unit sphere centered at the satellite, thus providing excellent directional sampling 
of the radiation intensity. 

The purpose of this paper is to analyze and describe the rigid body motion of 
this satellite, finally obtaining the angle that the counter axes make with the magnetic 
field lines as a function of time. This angle is then correlated with observed counting 
rates resulting in plots of counting rate versus angle that the detector makes with the 
magnetic field. 

/" 

The scope of this paper will be limited to obtaining these counting rate plots. 
The intensity as a function of angle may be determined by solving an integral equation 
containing these observed count rates and the angular response of the counter for 
various energies. This latter input awaits the completion of a new detailed calibration 
presently being undertaken by Baicy at the Ballistic Research Laboratories, Aberdeen, 
Maryland. 

The first chapter of this paper is of an introductory nature giving a brief re- 
sume and history of the radiation belts as well as some characteristics of the satellite 
pertinent to the analysis. This is followed by the determination of the total angular 
momentum vector configuration from antenna patterns. It is shown that body torques 
apparently are acting on the satellite to produce changes in the space-fixed components 
of the angular momentum vector. Verification of this determination is obtained by a 
confirmation of the anomalous temperature fluctuations and drag variations which are 
explained by this analysis. Finally the counter orientation is obtained at various times 
of interest and the corresponding counting rates are plotted as a function of angle be- 
tween the counter axis and the magnetic field. 
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A. THEORY OF PARTICLE TRAPPING IN MAGNETIC FIELDS 

An intensive study of the influence of the geomagnetic field on cosmic 
particles was undertaken by Stormer [Ref. i] and Alfven [Ref. 2] . While these works 
contain all the fundamental concepts necessary to understand the trapping of particles 
in the geomagnetic field near the earth, no obvious method for the injection of such 
particles was then apparent; hence Stormer and Alfven were primarily concerned with 
the motions of particles of extraterrestrial origin with high magnetic rigidities. Inter- 
est in magnetic trapping of charged particles was stimulated by attempts to confine 
high temperature plasmas with the hope of obtaining thermonuclear reactions. The 
magnetic mirror concept, although not by any means new (dating back to Poincare) , is 
presented very clearly by Spitzer [ Ref. 3] and would be worthwhile to review at this 
point. 


Consider a particle with mass m and charge q moving with velocity v in a uni- 
form magnetic field B. The basic equation of motion is 


d v — v 

m — = q (v x B) . 


( 1 ) 


If v_^is resolved into a component ^ perpendicular to B, and a component v parallel 
to B, it may be seen that the force is always perpendicular to both the velocity and 
magnetic lines of force. This results in the particle gyrating about a point on the line 
of force called the guiding center. Since the v |( component is unchanged, the guiding 
center will move along the line of force with velocity v (| . The resulting motion is a 
helical path about the magnetic field line. The radius of gyration a and the cyclotron 
frequency u> c are immediately found by equating the centrifugal force to the magnetic 
force exerted on the particle: 


or 


and 


m Vj_ 
a 


= q v. b, 


a = 


m v x 
q B ’ 


co=- = ^ 
c a m 


(2) 

(3) 

(4) 
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The magnetic moment of the gyrating particle is the current q co c /2 v times 
the fvoa encircled by the gyrating particle , ir a?. 


a 2 qco m v. 2 

~ 2 = 2B 


( 5 ) 


Consider the case of the spiraling particle moving into a region of greater B. 
Assume that the change B is small over a cross section of the particle's orbit such that 
9 B z /9z is constant in a cylindrical coordinate system with B along the z-axis , and is 
essentially equal to 9B/9z. From 


one obtains 


V . B = 0, 


i_ _9_ 
r 9r 


(r B r ) + 


^2 
9 z 


= 0 . 


( 6 ) 


(7) 


From the assumption on 9B z /9z, equation 7 may be integrated to give 


B r 


1 9 B 

2 9z 


( 8 ) 


Since B is along the z-axis, denote 9B/9z = B. Since B is increasing in density, 
the lines of force are converging and there must exist a radial component B r . This 
component through equation 1 results in a retarding force on the particle given by 

dv 

m 5T (9) 


From equation 8 with r = a. and equation 9 , 
dv 

m 5T = - 5 q 'l a ^ B = - < I0) 

Using this result with equation 5 and the assumption that kinetic energy is conserved, 
it may be shown that the magnetic moment p is conserved for slowly varying B ( See 
Appendix 1). This quantity is termed an "adiabatic invariant. " According to Northrop 
and Teller [ Ref. 4] the conservation of magnetic moment as well as other adiabatic 
invariants holds also for relativistic particles. 
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The conservation of magnetic moment has the immediate consequence that 
particles will be reflected away from a region of increasing magnetic field [Ref. 5] . 
This is seen from equation 5 if the velocity vector is assumed to make angle a with 
the line of force B. The component v. becomes v sin a and /u becomes 




mv* 

2B 


sin 2 a = constant. 


(ID 


As B increases, a also increases until sin a = 1. At this point, determined 
by the injection conditions B 0 and Q! 0 , the particle undergoes reflection and emerges 
from the region of denser B. In the geomagnetic field a particle injected with a 
magnetic moment such that mirroring will occur before the particle enters the denser 
atmosphere will continue to undergo reflections at mirror points in both hemispheres, 
thus becoming trapped on a line of force. Due to perturbations such as the perpendicu- 
lar gradient of B, centrifugal force due to particle motion along curved lines of force, 
and possible due to electric fields, the particles undergo a longitudinal drift. Assuming 
axial symmetry of the geomagnetic field, a trapped particle must remain on a figure of 
revolution described by the field line on which it is trapped. The nonaxial case is also 
discussed by Northrop and Teller. 


B. OBSERVATIONS OF TRAPPED RADIATION 


Prior to 1958 the only facilities available for high altitude measurement of 
cosmic rays were aircraft, rockets, and balloons. The intensities observed at 
altitudes available to these vehicles agreed quite well with expected intensities caused 
by cosmic ray primaries. In the studies of auroral activity using rocket-borne radi- 
ation detectors, Van Allen and others obtained evidence that the particles responsible 
for the aurora may be trapped on lines of force [ Ref. 6] . 

The launching of satellite 1958 Alpha, equipped with a State University of Iowa 
radiation experiment designed to detect cosmic primaries [ Ref. 7] , provided the first 
real evidence of the high intensity radiation zone beginning at 1000 kilometers above the 
earth at the magnetic equator. Actually in the case of 1958 Alpha, the intensity at times 
was so high that the counters were saturated, resulting in no observed counts [Ref. 8] . 

Despite this saturation effect, the observation that the radiation intensity in- 
creased rapidly in going from 700 km to 1000 km led to the conclusion that the high 
intensity radiation was constrained from reaching lower altitudes by the geomagnetic 
field; hence, the observed radiation must be charged particles trapped on magnetic 
lines of force. The correctness of this assumption was verified with the launchings of 
1958 Gamma and Epsilon or Explorers in and IV as they are commonly called [Ref. 9] . 

During the telemetry lifetime of Explorer IV the Argus Experiment, suggested 
by Christofilos[ Ref. 10] , was carried out. The experiment consisted of exploding three 
small yield atomic devices at different times at altitudes of around 400 km near the 
auroral latitudes. High energy electrons resulting from the burst and from beta decay 
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of the fission fragments became trapped on magnetic lines of force , thus forming an 
"Argus shell. " The configuration and intensity of this shell was then measured using 
Explorer IV [Ref. 11] and a series of Jason sounding rockets [Ref. 12] . 

Shortly after the Argus Experiment, the launching of several lunar probes 
made it evident that there existed a second radiation belt extending out to several earth 
radii [Refs. 13 through 16] . 

Other experiments, particularly those of Freden and White [Ref. 17] in which 
emulsion packs were flown and recovered in Thor nose cones , indicated that the ener- 
getic component of the inner belt was primarily due to protons. Also, it was found 
that the outer belt consisting of electrons is subject to temporal variations and is 
strongly connected with solar and auroral activity [Refs. 18, 19] while the inner belt 
remains relatively unvarying with time [ Ref. 20] . 

C. ANISOTROPY OF TRAPPED RADIATION 


For every point on a magnetic line of force there is some critical angle 
that the velocity of the particle must make with the field line if it is to be reflected 
before entering the dense atmosphere. Particles injected into the field with angles less 
than this critical value will soon be lost; hence, there should be a comparatively radi- 
ation-free cone defined by this critical angle. If a directional detector flown on a 
satellite is allowed to sweep through various angles with the field lines , strong vari- 
ations in counting rate should be observed as the counter moves in and out of this 
excluded zone. If the angle that the counter axis makes with the field line can be 
determined at various times , this may be correlated with the observed count rate at 
these times resulting in a determination of the directional distribution of radiation 
flux. Such a determination may be useful in determining the injection mechanism for 
the trapped particles [Refs. 21 through 23]. 


This angular distribution together with the integral invariant may provide a 
means for mapping the entire radiation belt using a theorem developed by Northrop, 
and Teller, that is, "In a steady state with no electric fields present, contours of 
constant B on an invariant surface are also contours of constant particle density [Ref. 


An effort was made to determine the angular distribution of trapped particles 
by Welch using the data from the Jason sounding rockets. The resulting plot of count 
rate versus angle very closely resembled the angular response of the detector; hence, 
it was concluded that the intensity function was indistinguishable from a delta function 
at 90 degrees to the magnetic field [Ref. 25] . These data were obtained at low altitude 
near a mirror point. 
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Van Allen reports that the Argus radiation was always observed to be disc-like 
with an angular thickness of the order of + 20 degrees or less [ Ref. 26] , and that the 
natural radiation appeared to be + 30 degrees [Ref. 27] . 

From an examination of the Channel 2 scintillation detector telemetry and the 
field strength recordings taken by the Huntsville Tracking Station, it was concluded that 
all the inf ormation necessary to obtain a detailed determination of angular distribution 
of radiation flux at various magnetic latitudes was available. The task of unscrambling 
this information and providing this determination was then undertaken by the Author. 

D. DESCRIPTION OF THE SCINTILLATION DETECTOR 

The scintillation detector on Explorer IV is described by Van Allen as con- 
sisting of a plastic scintillator disc, 0. 762 cm in diameter and 0. 178 cm thick. It is 
cemented on the face of an RCA 6199 photomultiplier tube [Ref. 28] . The axis of the 
detector was directed perpendicularly to the axis of symmetry of the satellite and 
through an open hole in the wall of the stainless steel payload shell, as shown in 
Figure 1. This aperture was covered by 0. 14 gm cm 1 2 of aluminum. The unidirectional 
geometrical factor* was 0. 040 cm 2 steradian, which increases for stopping power 
greater than 1. 6 gm cm 2 to an asymptotic value of 4. 2 cm 2 steradian. 

The collimation was such that the counter accepts the full value of incident flux 
of particles penetrating the aluminum window but not the steel shell for half angles up 
to 6 degrees. The acceptance falls linearly to zero at half angles of 19 degrees. 

The over-all response of the system may be represented by 


where t is the counter dead time of 91 ju sec, r] is observed counting rate, and N is true 
counting rate. 

This detector responds to electrons greater than 680 kev, protons greater than 
10 Mev, and X-rays greater than 300 kev but with low efficiency [Ref. 29] . 

The output from this detector is scaled by a factor of 2048 and then fed to the 
Channel 2 bistable subcarrier oscillator. At a signal from the scalar, the subcarrier 
changes from the high frequency state ( 602 cps) to the low frequency state ( 518 cps) . 
The output from the discriminator at the receiving station is seen as a square wave 


1. This factor is defined by Van Allen as the true counting rate divided by the uni- 

directional intensity (particles cm -2 sec -1 steradian -1 ) and the efficiency of the 
detector for the particles in question. 
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& HIGH POWER TRANSMITTER 

FIGURE i. EXPLORER IV PAYLOAD CONFIGURATION SHOWING THE PLACEMENT 
OF THE DIRECTIONAL SCINTILLATION DETECTORS. THE CHANNEL 
2 COUNTING RATES ARE THE OUTPUTS FROM DETECTOR A 




with a period representing the time required to collect 2048 counts. Figure 2 is a 
sample of telemetry which clearly shows the effect of the satellite body motions on the 
Channel 2 counting rate. 

E. GENERAL DESCRIPTION OF THE RIGID BODY MOTION OF SATELLITE 
1958 EPSILON 

The payload configuration of 1958 Epsilon is shown in Figure 3. For the 
purpose of discussing the rigid body motion, the satellite may be approximated by a 
long cylinder with a longitudinal to transverse moment of inertia ratio of 1 to 84. 

Since the upper stages of the Jupiter C launching vehicle were unguided, the 
payload and the upper stage cluster were spun before injection to assure stability dur- 
ing the coast phase of the trajectory. This resulted in the payload having an initial 
spin of 24 7 r rad/sec about the axis of symmetry which is the axis of least moment 
of inertia. Since the rotation is about the least moment of inertia, the satellite is in a 
maximum rotational energy configuration which is unstable for a semi-rigid body. 

In the case of the earlier Explorers (I and HI) , a transition was observed from 
the axial spin at injection to a tumbling motion about a transverse axis [Ref. 30] . A 
detailed analysis of this transition is given by Wells [Ref. 31] for the case of Explorer 
I which had flexible wire turnstile antennas. This transition may be understood by 
considering a slight malalignment of the angular momentum and angular velocity 
vectors at injection. The axis of symmetry would then perform force-free precession 
about the angular momentum vector with precession rate <£. This precession rate may 
be represented by means of a Poinsot construction as outlined by Goldstein [ Ref. 32] , 
or by the solution of Lagrange T s equations expressed in terms of the Euler angles <p , 

0, and ip. This latter treatment is given by Slater and Frank in which the z-axis is 
taken as the direction of the total angular momentum vector, 0 is the angle the figure 
axis of symmetry makes with the z-axis, <p is the angle the line of nodes makes with 
the space-fixed x-axis, and ip is the angle the body-fixed x-axis makes with the line of 
nodes. Let the transverse moments of inertia be equal and designated by Ij and let the 
longitudinal angular momentum be I 3 o> 3 . From equation 3. 4 in Slater and Frank [Ref. 
33] 

h 

<p cos 0 = — ct> 3 . ( 13) 

h 

This derivation is shown in Appendix 2. From the definition of the coordinate system, 
the total angular momentum is related to I 3 co 3 by 

L cos 0 = I 3 u >3 . ( 14) 
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FIGURE 3. EXPLORER IV EXTERNAL CONFIGURATION AND DYNAMICAL 
PROPERTIES 



The precession rate <p is then found to be 


0 


L_ 

Ii 


(15) 


Initially all the angular momentum is manifested in rotation about I 3 ; hence, 

<p = — oo 3 (initial) . (16) 

Ii 


With any small precession, a body-fixed mass. element will experience oscil- 
latory driving forces because of the precession rate <p and body-fixed roll rate ip. 

Unless the body is completely rigid, this mass element will be driven in forced 
oscillation resulting in a dissipation of rotational energy through frictional losses. 

Since no external torques have acted on the body, L must be a constant of the motion 
while rotational energy is lost. Therefore, the body must rotate about a greater 
moment of inertia as u> 3 decreases. Since from equation 14 <fi is constant as long as 
L is conserved, it is seen from equation 13 that cos 6 decreases with co 3 . Con- 
sequently, the opening angle or the angle between the axis of symmetry and L increases 
with' decreasing energy. Eventually this angle approaches 90 degrees, the roll rate 
goes to zero, and the driving forces vanish. The body is then in the stable or minimum 
energy configuration, and is tumbling, propeller-like, about the initial angular mo- 
mentum vector which will also be referred to as the tumble axis. 

The rate at which this transition takes place may be calculated if the dissipation 
function is known. In the case of Explorer I, the whip antennas were by far the largest 
contributor to the energy dissipation and were successfully treated as damped pendu- 
lums. However, in Explorer IV, the dissipation mechanism is not so obvious since ef- 
forts were made to make this payload as completely rigid as possible. Nevertheless, 
after several days of orbiting, evidence of tumbling was observed as the antenna power 
pattern nulls along the nose and tail [ Ref. 34] produced periodic fades in the recorded 
field strength at 3. 5-second intervals. Since the moment of inertia was 1 to 84 and 
initial roll rate was 24 n rad/ sec , this fading period corresponded to half the preces- 
sional period. 

Since the antenna pattern was virtually symmetrical about the polar axis, it was 
not possible to obtain roll rate from radio observations. However, roll rate could be 
ascertained from the disc-like character of the trapped corpuscular radiation and the 
counting rate modulations of the directional counter. Such a study was performed by 
Fields [Ref. 35] resulting in the roll period history shown in Figure 4. From this 
result and equation 13, it may be seen that the opening angle had become 80 degrees 
after five days of orbiting. Note also that the roll never goes to zero; in fact, it is not 
even a monotonic decreasing function. The anomalous fluctuations on August 28 and 29 
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JULY 1958 


AUG. 1958 


SEPT. 1958 


FIGURE 4, 


OBSERVED ROLL PERIOD DURING THE TELEMETRY LIFETIME 
OF EXPLORER IV 



may be due to external torques which somehow have come into play. Even with a 
slight residual roll rate, the roll angular momentum can be neglected compared to the 
tumble angular momentum, and for all practical purposes the opening angle can be 
considered to be 90 degrees after August 1. 

To determine the counter orientation, the first task is to specify the tumble 
axis in some inertial coordinate system. At injection, the angular momentum vector 
is directed along the velocity vector whose components are known in an inertial refer- 
ence frame; hence, the initial tumble axis is given. Under the assumption that no 
external torques act on the satellite, the orientation of tumble axis would remain fixed 
in space. However, evidence to the contrary is reported by Snoddy through measure- 
ment of the temperature of the satellite interior [Refs. 36, 37] . Figure 5 is a history 
of the interior temperature obtained from the high and low frequency states of the 
Channel 3 subcarrier oscillator. The temperature dependence of this subcarrier oscil- 
lator frequency was determined by preflight calibration by Ludwig at the State Uni- 
versity of Iowa. The initial drop in temperature may be attributed to the transition 
from axial spin, with the sun at almost a right angle to the symmetry axis, to the state 
of tumble where the insolated area time averaged over a tumble cycle is only 63. 7 per 
cent of the longitudinal cross-sectional area. This decrease in insolated area results 
in less heat input and a lower equilibrium temperature. The fluctuations in tempera- 
ture after the initial drop are difficult to explain except in terms of varying insolated 
area. Knowledge of the thermal properties of the satellite and of the position of the 
orbit relative to the earth and sun allow calculation of the expected temperatures for 
various percentages of insolated areas. The results of these calculations performed 
by Snoddy are shown in Figure 6. The maximum and minimum measured temperatures 
fall almost within the limits of the expected values for insolated areas ranging from 1 
to . 65 of the longitudinal cross section. This is strong evidence that the orientation of 
the plane of tumble is changing relative to the sun. The motion of the sun can account 
for a change of approximately one degree per day, but the area change needed to ex- 
plain the temperature data demands an orientation change of 90 degrees in about 9 
days, or an average change on the order of 10 degrees per day. 

Throughout the telemetry lifetime of the satellite, the tumble period was ob- 
served to be very close to 3. 5 seconds. Therefore, it appears that angular momentum 
is conserved in magnitude, but that through some mechanism, non-dissipative external 
torques are acting on the satellite to produce a change in the orientation of the angular 
momentum vector. 
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FIGURE 5. EXPLORER IV TEMPERATURES MEASURED DURING TRANSMITTER 
LIFETIME FROM THE CHANNEL 3 SUBCARRIER FREQUENCY 
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FIGURE 6. CALCULATED TEMPERATURES FOR EXPLORER IV FOR VARIOUS 
FRACTIONS OF THE LONGITUDINAL AREA PRESENTED TO THE 





SECTION II. DETERMINATION OF THE TUMBLE AXIS 


A. DISCUSSION OF THE TUMBLE AXIS DETERMINATION 

As previously pointed out, the orientation of the tumble axis must first be 
determined in order to obtain the counter orientation. This becomes somewhat pro- 
blematical since the satellite was not provided with any onboard sensors for this pur- 
pose. However, examination of the satellite antenna radiation patteims suggests a 
method whereby the look angle, or the angle that the line of sight from the tracking 
station to the satellite makes with tumble axis, may be inferred from the recorded 
signal strength variations [Ref. 38] . Figure 7 shows these radiation patterns for the 
108 mb/sec and 108. 03 mc/sec transmitters. It may be seen that the 108 mc/sec 
antenna, being an asymmetrical dipole, has a minor lobe on either side of the nose of 
the satellite whereas the 108. 03 mc/sec antenna exhibits a characteristic dipole pattern. 
In either case the radiation is linearly polarized in a plane containing the symmetry 
axis. 


If the satellite is viewed along the tumble axis with a linearly polarized re- 
ceiving antenna, fading will be observed at twice the tumble frequency because of the 
rotating plane of polarization of the emitted radiation. The plane of polarization is 
also rotated by the Faraday effect resulting from magneto-ionic splitting [Refs. 39, 40] 
as the signal traverses the ionosphere. The rotation rate of the latter effect, which is 
due to the changing integrated electron density along the path of propagation as the 
satellite moves relative to the tracking station, is' slow compared to the tumble rate 
for frequencies of 108 mc/sec. Therefore, it does not materially change the observed 
patterns of signal strength variations. Thus, for the case where the look angle is 
small, that is, when the satellite is viewed almost along the tumble axis, the signal 
strength variations are essentially the same for both transmitters. 

If the satellite is viewed in the plane of tumble, fades will occur in the 108. 03 
mc/sec signal strength patterns at twice the tumble frequency due to the nulls in the 
antenna radiation pattern. These nulls are indistinguishable from those caused by 
polarization fading in the case of small it'; consequently, the 108. 03 mc/sec signal 
strength variations do not provide usable information if a linearly polarized tracking 
antenna is used. However, when viewing in the tumble plane, the small lobes on 
either side of the nose are apparent in the 108 mc/sec signal strength patterns. This 
distinguishing feature then enables one to determine whether the look angle \I> is large 
or small. As it turned out, most of the signal strength recordings taken by the Hunts- 
ville Tracking Station contained the small lobes in the 108 mc/sec pattern. On two 
days, August 12 and 18, passes were observed in which the small lobes were absent 
for an interval of time during the pass. A simulation of the observed signal strength 
variations was performed on an analogue computer in which it was shown that the small 
lobes would be visible in' the signal strength variations unless ^ was less than about 
20 degrees. 
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FIGURE 7. EXPLORER IV ANTENNA RADIATION PATTERNS 



The line of sight vector was computed at the times where the small lobes were 
absent on the August 12 and 18 passes. This computation was done by subtracting the 
vector from the center of the earth to the tracking station from the radius vector to 
the satellite. The former is known from the geodetic coordinates of the station. The 
latter is found from the geodetic position of the satellite which is listed in one-minute 
intervals in the Smithsonian Astrophysical Observatory emphemeris [Ref. 41] . The 
resultant line of sight vector was then transformed into an earth-centered, right- 
handed, space orientation-fixed coordinate system in which the x-axis is in the direc- 
tion of the vernal equinox and the z-axis is along the north celestial pole. To facilitate 
the plotting of this vector, which was normalized and designated as U , it is desirable 
to express it in terms of two position angles: right ascension and declination. Since 
the satellite is tumbling essentially in a plane, one should expect to observe the same 
signal strength variations on either side of the plane. Hence, the same pattern should 
be observed when looking along U as when looking along -U. Therefore, the U was 
plotted as well as -U for each pass. Figure 8 is such a plot of the sight vector for 
two passes on August 12. The single lines indicate the orientations of the sight vector 
when the small lobes were visible on the signal strength recordings - hence ^ > 20 
degrees. The double lines indicate the orientations of the sight vector when these 
lobes were absent or when it < 20 degrees. The dashed line around the region where 
iff was observed to be less than 20 degrees represents the intersection of a cone whose 
half angle is 20 degrees and the celestial sphere. This fit is determined graphically 
and the axis of the cone is assumed to be parallel to the tumble axis. At this point, 
there is no way of determining the sense of tumble; hence, the angular momentum vec- 
tor of the satellite may either be parallel or anti-parallel to the cone axis. Figure 9 
is a similar determination of the tumble axis obtained from two passes on August 18. 

Unfortunately, many passes recorded by the Huntsville Tracking Station prior 
to August 18 did not contain useable signal strength recordings. Also, the condition 
that iff be less than 20 degrees to provide useful orientation data was not frequently 
met. It became apparent that if a history of the tumble axis orientation over the satel- 
lite transmitter lifetime was to be determined, an alternate method must be found. 

On August 18, the Huntsville Tracking Station altered' their receiving antenna system. 
Instead of using a single dipole array, an attempt to obtain polarization diversity by 
adding another dipole array with crossed polarization was made. The two antenna 
arrays were mounted on separate pedestals, being separated by about 50 meters. 

Their outputs were added a "T" junction arid then detected by a single Hallamore re- 
ceiver. As the satellite moved overhead, the path length to the two antennas changed, 
thus altering the phase angle between the two antenna outputs. This produced an 
interferometer effect that altered the effective polarization of the receiving antenna. 

For example, if the two output voltages are in phase, the array is effectively plane 
polarized at 45 degrees to each element. As the phase angle becomes 90 degrees, the 
array becomes circularly polarized. The location of the antennas was such that the 
polarization change was usually slow compared to the tumble rate so that several 
tumble cycles could be observed at more or less constant polarization. 
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FIGURE 8. RIGHT ASCENSION AND DECLINATION OF SIGHT VECTOR FOR 
PASSES DURING AUGUST 12, 1958 






This technique allowed the use of 108. 03 mc/sec signal strength patterns. 

The times when the antenna array was circularly polarized could be obtained by ob- 
serving when the fades were less severe. At these times, the amplitude of the fades 
allowed an estimate of the look angle. Lack of calibration on the recorded signal 
strength records precluded an accurate measurement of the look angle; but a crude 
estimate could be made on the basis of comparison of the fading amplitude observed 
during the circularly polarized portion with the more severe fades observed when 
the antenna was linearly polarized. If the signal strength exhibited periods where 
little or no fading was evident, such as shown in Figure 10, the look angle was 
estimated to be less than 30 degrees. If the fading was discernible, but still much 
less than the severe fades as in Figure 11, ^ was estimated to be between 30 degrees 
and 45 degrees. If the observed signal was a mixture of antenna radiation pattern 
fades and polarization fades as shown in Figure 12, ^ was estimated to be from 45 
degrees to 60 degrees. If the pattern was predominately antenna radiation pattern 
fades, lessened only slightly as the polarization changes, was estimated to be 60 
degrees to 75 degrees. Figure 13 is an example of this pattern. If was between 75 
degrees and 90 degrees, severe fades due to the radiation nulls along the nose and 
tail were observed regardless of the polarization of the receiving antenna. This is 
shown in Figure 14. 

For days when a sufficient number of passes were available, the right ascension 
and declination of the 1 sight vector at various times during the passes were plotted in 
the manner described previously. The type of signal strength variation noted is indi- 
cated by shading. A set of transparent overlays were prepared showing the inter- 
sections on the celestial sphere of a family of cones of various half angles about an 
axis having a particular declination. These overlays were used to determine graphi- 
cally the position of the tumble axis that best satisfied the observed signal strength 
patterns. Figures 15 to 20 show this determination for various days in August and 
September. The right ascension of the tumble axis can be fairly well determined by 
making use of the 75 degrees to 90 degrees patterns in the signal strength which 
appear to be fairly reliable. The declination is probably determined with less cer- 
tainty since it is generally near the equator and not well bracketed by the 75 degrees to 
90 degrees patterns. 

An effort was made to obtain more information on the orientation of the tumble 
axis by examining data from other tracking stations. Unfortunately, little usable data 
was found since many stations did not record signal strength or did not use antenna 
systems that facilitated analysis of the signal strength variations in terms of look 
angle. Some useful data was obtained from the station operated by the Stanford Re^ 
search Institute located at Johnston Island. This station used a steerable, circularly 
polarized, crossed Yagi antenna; however, they did record signal strength. It was 
noted in the 108. 03 mc/sec telemetry records taken by this station that were intervals 
where the count rate data was uninterrupted by fades. At other times, the telemetry 
was interrupted every 3. 5 second by fading that was presumably caused by the signal 
strength dropping below the receiver sensitivity. When ^ is small the signal strength 
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FIGURE 11. TELEMETRY SAMPLE ILLUSTRATING THE SIGNAL STRENGTH 
VARIATIONS FOR THE CASE OF 30° < * < 45° 
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FIGURE 16. RIGHT ASCENSION AND DECLINATION OF SIGHT VECTOR FOR 
PASSES DURING AUGUST 27, 1958 
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FIGURE 17. RIGHT ASCENSION AND DECLINATION OF SIGHT VECTOR FOR 
PASSES DURING AUGUST 30, 1958 
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variations are small, and since the antenna is not affected by polarization of the 
incident radiation, one would expect uninterrupted data. When ’F is large, fades would 
be expected because of the observance of the radiation nulls along the satellite nose and 
tail. It is not known what values of 'k will just begin to produce fades since this de- 
pends on the receiver sensitivity, the distance from the satellite to the station, the 
power radiated by the satellite and any attenuation suffered along the path of propa- 
gation. The orientation of the sight vector was plotted at the times the signal was 
uninterrupted resulting in Figures 21 and 22. The dashed line indicates a cone of 45 
degrees half-angle which appears to define the region on the celestial sphere in which 
the signal strength does not drop below the receiver sensitivity. The axis of this cone 
is again assumed to be parallel to the tumble axis. 

The methods used for determining the tumble axis are admittedly crude, since 
the look angles could only be estimated. Graphical fitting was used primarily for 
expediency. While a more elegant numerical method could have been attempted, it 
was not felt justifiable in light of the crudeness of the data. Furthermore, it was 
necessary to use several passes during the same day to define the tumble axis suffi- 
ciently. Since the tumble axis may have shifted by as much as 10 degrees during this 
time, the best that can be hoped for is some mean location that may be considered 
effective for that day. 

To obtain an analytical description of the excursions of the tumble axis, the 
right ascensions and declinations obtained by this analysis were plotted as a function 
of days after launch. A polynomial description using two 9th-degree polynomials was 
attempted. Figure 23 shows the result of attempting to fit polynominals to the right 
ascension and declination of the tumble axis. No data is available from launch until 
August 9; hence, the polynomials are not defined except from August 9 to September 11. 
The polynomial description for declination from September 3 to September 9 does not 
appear to fit the observed data very well, but because of the previously mentioned un- 
certainties in determining declinations near the equator this discrepancy will be dis- 
regarded. 

B. CONFIRMATION OF THE TUMBLE AXIS DETERMINATION 

Having obtained polynomials that describe the motion of the tumble axis, the 
angle made by the position vector of the sun with the tumble axis may readily be found. 
From this angle the average insolated area during a tumble cycle may be calculated 
Then using the curves of predicted interior temperature for various insolated areas 
prepared by Snoddy (Fig. 6) , the temperature may be found as a function of time. 

This result is plotted over the observed temperature in Figure 24. The dotted line in 
the plot represents the temperature predicted if the tumble axis had retained its initial 
spatial orientation. There is a discrepancy in the temperature ^obtained from the high 
frequency state of the subcarrier oscillator as compared to that obtained from the low 
frequency state which is probably due to some frequency drift which occurred after the 
calibration was made. This discrepancy appears to be uniform throughout the 
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FIGURE 21. RIGHT ASCENSION AND DECLINATION OF THE SIGHT VECTOR FOR 
AUGUST 9, 1958 AT TIMES WHEN FADING WAS NOT OBSERVED IN 
THE TELEMETRY RECEIVED BY THE JOHNSTON ISLAND STATION 
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FIGURE 22. RIGHT ASCENSION AND DECLINATION OF SIGHT VECTOR AT 

TIMES FADING WAS ABSENT FROM JOHNSTON ISLAND TELEMETRY 
ON AUGUST 11, 1958. 
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transmitter lifetime. Although the actual temperature may not be known exactly, the 
observed temperature fluctuations appear to be consistent with those predicted from 
the determination of the tumble axis. 

The change in orientation of the plane of tumble should also result in drag 
variations because of a change in average frontal area presented to the atmosphere at 
perigee. Such fluctuations were observed by'Lautman at the Smithsonian Astrophysi- 
cal Observatory [Refs. 42, 43] . Figure 25 shows the period loss per revolution due 
to drag forces. Using the polynomials for the tumble axis description, the angle that 
the tumble axis makes with the perigee velocity vector was found and the resultant 
average drag area was calculated. Figure 26 shows the computed average frontal 
area as a function of time. It may be seen that, in general, the loss of period is 
greatest when the drag area is large. Of course there are other factors that influence 
drag losses, such as a variation of atmospheric density which appears to be related 
to solar activity [ Ref. 44] . 

In a study of the effect of the tumble axis orientation relative to magnetic field 
on counting rate observed by the directional radiation detectors on Explorer IV, 

Shelton points out that certain counting rate patterns should be characteristic of 
certain values of the dot product of the tumble axis L and the magnetic field direction 
B [Refs. 45, 46] . The polynomial for L was incorporated in the program [Ref. 47] 
for computing the satellite position and the local geomagnetic field components using 
the Finch and Leaton [Ref. 48] 48 term expansion. The result of this program was a 
listing of the dot product of L and B at one-minute intervals from August 9 through 
the transmitter lifetime. Using these results, Fields was able to match the observed 
counting rate patterns with those predicted by the computed L and B on several occa- 
sions [ Ref. 49] . 

It is difficult to say how accurate the previous tumble axis determination actu- 
ally is, although the success in explaining the observed temperature and the drag 
variations is encouraging. It is felt that the determination is generally good to about 
+ 10 degrees in angular position with possibly somewhat higher uncertainty in declina- 
tion during the times when the tumble axis lies near the equator. It will be shown 
later that uncertainties of this magnitude do not seriously detract from the accuracy of 
obtaining count rate as a function of the angle the counter axis makes with the magnetic 
field lines. 
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SECTION in. DETERMINING COUNTER ORIENTATION 


It was shown by Shelton that the angle the counter axis makes with the magnetic 
field may be conveniently expressed in an appropriate choice of coordinate systems 
[ Ref. 50] . Let A be a unit vector along the satellite symmetry axis and P be a unit 
vector along the counter axis. In a double-primed coordinate system, let A" be along ! 

the y" axis. Since from the geometry of the satellite (Fig. 1) P is perpendicular to 
A, P will lie in the x" - z" plane making angle (ip t + 6) with z". The counter axis ? 

may then be expressed by 



where ip is the roll rate and 6 is the roll phase angle or the angle P makes y" at 
t = 0. 

Let B be a unit vector along the local geomagnetic field at t = 0. Consider a 
primed coordinate system defined such that x' is along the tumble axis L and B lies 
in the x' - y' plane. Since A is perpendicular to L, the symmetry axis is given by 


A' = 


cos (<p t + y) 
sin (<p t + y) , 


(18) 


where (p is the tumble rate and y is the tumble phase angle or the angle A ' makes with 
y' at t = 0. 

The counter axis P may be found by 


1 0 0 \ / sin (ip t + 6) 

P ' = ( 0 cos (<pt + y) - sin ( <p t + y) I I 0 

0 sin (<pt + y) cos ( <p t + y) 1 \ cos (ip t ’+ d) 


(19) 


Let B' have components bj, b 2 , b 3 along x' , y', and z', respectively. The 
angle A between the counter axis and the magnetic field is found immediately by 


cos A = P' ■ B' = b 1 sin (^ t + 6) - b 2 sin (0 t + y) cos (ip t + 5) 
+ b 3 cos ( <p t + y) cos (ip t + 5) . 


( 20 ) 
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From the equation 20, an idea of the expected periodicity of the counting rate 
variations may be obtained. The choice of reference time or t = 0 is completely 
arbitrary. At t = 0, b 3 = 0 f rom the definition of the primed coordinate system. Then 
bj becomes L • B and is Vl - b A 2 . If the tumble axis is along B, the b 2 term van- 
ishes and A depends only on ip t. This results in the counting rate being modulated at 
the roll frequency. However, if L is perpendicular to B, the b A term vanishes and the 
counting rate will be modulated by the $ term when cos ( ip t + 6) « i. Therefore, 
periods of tumble modulations will be observed until cos (ip t + 6) becomes small. At 
these times P remains nearly perpendicular to B and the counting rate is more or less 
constant. Such modulations have been observed and pointed out by Van Allen [Ref. 51] 


From this expression the angle A may be found at various times during a parti- 
cular passJLf <p , ip , y, 6 , b l9 b 2 , and b 3 are known. The b*, t^, and b 3 may be found 
from the L and B by 



( 21 ) 


where B x , By, B z are the normalized components of the field in a vernal equinox co- 
ordinate system and K = L * B. The Smithsonian ephemeris lists the satellite position 
and the local geomagnetic field components in one-minute intervals throughout the 
transmitter lifetime. These magnetic field components are expressed in gauss and 
resolved in terms of a local vertical system in which B D is antiparallel to the local 
vertical, B E is perpendicular to the local vertical and extends to the east, and B N is 
mutually perpendicular to B D and B E and is directed to the north. These components 
may be transformed into the vernal equinox system by 



( 22 ) 
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where R x , Ry, R z are the direction cosines of the radius .vector from the center of 
the earth to the satellite in the vernal equinox coordinate system. The quantities Bn, 
Bg, and Bp are the normalized components of the magnetic field as given in the 
Smithsonian ephemeris. 

The tumble rate and phase may be found by making use of the nose and tail 
radio frequency radiation pattern nulls observed in the signal strength recordings. As 
discussed in Chapter It, the nose nulls are identifiable by the small lobes that appear 
on either side of the null. The observation of the null in the radio frequency radiation 
pattern along the nose or tail of the satellite occurs when the axis of symmetry makes 
the closest approach to the sight vector U. The axis of symmetry at the observance 
of a null is given by 


t* ( L x U) x L 

A = + — 

| L x U | 

where the minus sign is taken for a nose null and a plus sign for a tail null. 


(23) 


The angle through which the satellite has tumbled between the time of the first 
null and the time of the i^ 1 null is 



- (p Q = 2ii7T + arccos (Aj 



(24) 


where t^ is the time of the i^ 1 null after the reference time and where n is the number 
of revolutions the satellite has made between the first and the i^h null. The angle <p 0 
is the angular displacement of Aj at t = 0. The sign of the arccosine term depends 
on the sense of the tumble and the angular component of the U. 


The times at which nose or tail nulls were observed were then found during a 
particular pass, the A^ at each null was computed, and the angle through which the 
satellite had tumbled since the first null was found. If the tumble rate is assumed 
constant during a pass, the angular displacement must be a linear function of time. 
Imposing this requirement, the angular displacement was fit to a first degree poly- 
nomial by the method of least squares. The constant term is therefore the initial 
angular displacement <p Q and the first power coefficient is the angular rate <p. 

The accuracy of this determination should be emphasized. The primary source 
of error in this measurement is the time of the observance of the null. Since a time 
signal from WWV was recorded on the strip chart containing the telemetry, these 
times are probably accurate to . 1 second. The arccosine term in equation 24 is 
generally less than . 001 revolution per tumble cycle. Therefore, the uncertainties in 
L are not significant. The error in <p resulting from an error in t is 
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• A t 
A $ =$ — 


(25) 


so that the error is reduced by measuring over a long interval t. In general, the time 
available for observing nose and tail nulls is on the order of 5 minutes. An uncertainty 
of . 2 seconds in determining the interval t between the first and last null corresponds 
to i part in 1500. Since the tumble rate 0 is approximately a radian per second, the 
maximum error in 0 is a milliradian per second. Since least square fitting was 
employed, the error in 0 is probably on the order of + . 0001 rad/sec. Figure 27 
shows the values of 0 obtained from August 27 through September 1. It may be seen 
that the tumble rate is slowly decaying, having lost approximately . 008 rad/sec in 5 
days. This loss can be attributed to damping by both the atmosphere and the geo- 
magnetic field. However, the decay rate corresponds to a loss of 10 -5 rad/sec during 
a pass and does not invalidate the assumption that 0 is constant during a pass. 

This same accuracy is also retained in _Jhe phase angle 0 O . From this phase 
angle, the components of the symmetry axis, A, at t = 0 may be found by 



The determination of 0 and y thus far is completely independent of any as- 
sumption on the distribution of trapped particle flux. There is, however, still the 
ambiguity in the sense of tumble. Examination of the counting rate data reveals two 
important observations. During times that L is nearly perpendicular to B, sharp 
peaks in counting rate are observed at 3. 5-second intervals. Also, the peaks in 
counting rate are immediately preceded by a nose or tail null. If a particle is trapped 
on a line of force, the acute angle a that its velocity vector makes with a line of force 
as it approaches the satellite should be independent of the direction in which the guiding 
center moves along the line of force. Therefore, the distribution of trapped particle 
flux should be symmetrical about a plane perpendicular to the magnetic field. The 
angle Ap that the counter axis makes with the magnetic field when a sharp counting peak 
is observed may be found by neglecting bj and b 3 in equation 20, or 
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-cos Ap - sin ( cp X + y) cos {ip t + d\ 


(28) 


In the interval of interest, <p » ip , and cos {ip t-: + 6) may be considered constant dur- 
ing a cycle of <£. From the argument that the flux distribution is symmetrical about a 
plane perpendicular to B, a counting rate peak must be observed at + cos Ap as well 
as at - cos A^. In general, four such peaks should be observed during a tumble cycle 
unless Ap is 0 or 7r/2. Since particles whose velocity vectors lie along lines of force 
are not trapped, and since only 2 counting rate peaks are observed during a tumble 
cycle, it is concluded that the peak intensity is observed when A = 7r/2^. This results 
in counting rate peaks occurring when the symmetry axis is along the magnetic field 
line for the case where B is nearly perpendicular to L. This result, together with 
the observation that a nose or tail fade occurs just prior to a counting rate peak, 
allows the direction of the tumble to be determined. The geometry of the passes dur- 
ing the period of interest (August 27 through September 1) was such that the satellite 
rose above the radio horizon just after crossing the equator. The tumble axis was in- 
deed almost perpendicular to the line of sight and to the magnetic field, and was almost 
along the line of equinoxes. Figure 28 shows the relationship of the line of sight and B. 
Since a nose or tail null was observed as the symmetry axis passed the line of sight 
just prior to the observance of a peak counting rate as the symmetry axis passes U 5 
it is concluded that the angular momentum vector must be in the direction of the vernal 
equinox. 

The task of determining the roll rate and phase remains. Since it is not 
possible to infer roll from the signal strength, it becomes necessary to make use of 
the counting rate data. In the first case where b t « t> 2 , 

| cos A | — [ b 2 sin {<£ t+ y) - b 3 cos {<p t + y) ] cos {ip t + 6). (29) 

During the period of continuous high counting rate extending over several tumble cycles, 
it is presumed that A remains nearly ir/2 independent of <fi. This may happen only if 
cos {ip t + 6) becomes small or ip t + <5 ~ (2n-i) 7r/2 % - The times of such occurrences 
during a pass are noted, and the angular displacement of ip t + 6 is taken to be (2n-l) 
7r/2 at the times of the midpoints of such occurrences. Again employing the assump- 
tion of uniform rotation during the pass, the angular displacements are fit to a first 
degree polynomial by the method of least squares. In this manner the rate ip and the 
phase may be readily found. 


1 This does not necessarily imply that the intensity of particle flux is maximum at 
90° to the field. Since the counter has a finite resolution, a peak counting rate at 
A = 90° may be observed if the peak in particle flux is sufficiently close to 90° to be 
within the resolution of the detector. 
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FIGURE 28. GEOMETRY OF A TYPICAL PASS AT ACQUISITION DURING THE 
PERIOD FROM AUGUST 27 TO SEPTEMBER 1 



If bj is not small compared to be b 2 , ip and <5 may be determined by observing 
times of peak counting rates and equating cos X to 0 in equation 20 at these times. The 
resulting equation may be solved for ip t + <5 giving 

b 2 sin (0 t + y) - b 3 cos (<p t + y) 

I tan ( ip t + 6)1 = — (30) 

i b t 

From several counting rate peaks, various values of ip t + 6 may be found at various 
times. Again, # the least square fit of this data to a first degree polynomial may be 
used to obtain ip and 6. 

Again there is an ambiguity in the sign of ip. As long as the b* term remains 
small, cos X depends only on cos ( ip t + 6) and the sign of (ip t + 6) is immaterial. 

This is generally the case in the passes analyzed. The calculations were performed 
assuming both positive and negative roll directions at times when bj was the largest. It 
appeared that the positive value for roll yielded a better fit to the radiation data, but 
bj never became large enough to determine this definitely. 
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SECTION IV. RESULTS AND CONCLUSIONS 


A. ANALYSIS OF OBSERVED RADIATION 

Having developed a procedure in the previous chapter whereby the counter 
orientation may be determined as a function of time, particular passes during the 
period from August 27 to September 1 were selected for analysis. This period was 
chosen because it included the first two Argus events and the records were of ex- 
ceptional quality. Unfortunately , the Channel 2 subcarrier oscillator ceased opera- 
tion on September 2, shortly before the third Argus event; hence, no directional counter 
data is available after this date. 

In selecting the passes for this analysis, several requirements had to be con- 
sidered. First and most important, the radiation counting rate had to be high enough 
so that the accumulation time of the counter was short compared to a tumble cycle. 

In general, two passes each day met this requirement. Second, the pass must have 
readable time scales. Third, the pass must have readable radio frequency signal 
strength recordings in order to determine the tumble rate and phase. Eight passes 
were selected during the time from August 27 to September 1 that met these require- 
ments. They are passes 414, 415, 416, 428, 453, 454, 466, and 479.^ 

The Channel 2 counting rate was read using a Universal Telereader in con- 
junction with an IBM card punch. The distance on the telemetry record from a refer- 
ence point to each switch of the subcarrier oscillator was read and punched into an 
IBM card. In a similar manner the distance between two known times was recorded 
on an IBM card. These cards were used in a computer program that converted the 
location of a switch of the subcarrier oscillator frequency into time after the refer- 
ence point, determined the time interval between successive switches, and divided this 
time into the counter scale factor. The resulting count rate is taken to be effective at 
the midpoint of the interval between the successive switches. This time and the 
corresponding counting rate were then recorded on magnetic tape. To prevent cumula- 
tive time errors, a new reference time is taken every 15 seconds during the pass from 
the WWV time scales. Further error analyses follow in part B of this section. 

The times of occurrences of the signal strength fades because of the nose and 
tail nulls were determined and the (p , ijj , y, and <5 were calculated as outlined in 
Section III. These inputs along with the reference time and the components of L were 


1. The pass numbering convention used here takes the first ascending node as the be- 
ginning of Pass 1. Passes are then numbered successively at each subsequent as- 
cending node. 
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then fed into a computer program which contained the orbital elements and the coef- 
ficients for the magnetic field expansion. The magnetic tape containing times and the 
corresponding counting rates was then read by the computer and the angle X for each 
time was computed and listed with the corresponding counting rate. 

Figures 29 through 31 show the observed counting rate as a function of time 
together with the computer A. as a function of time. Note the counting rate peaks that 
are due to tumble modulations appear to coincide quite well with the X = 90 degrees 
points. It should be remembered that the tumble rate and its phase determination is 
completely independent of the radiation counting rates except for the geometric argu- 
ment for the direction of rotation. 

Figures 32 through 37 are plots of the observed counting rates as functions of 
X. Some scatter of points is observed and data smoothing techniques must be employed 
to find the best fit to this data. The degree of uncertainty associated with these points 
will be discussed in the following section. The smoothed data may then be corrected for 
counter dead time according to equation 12 to obtain the true count rate as a function 
of X. 


Examination of the counting rate versus X plots reveals a plateau of almost 
uniform low counting rates from X = 0 degree to about 50 degrees. This is probably due 
to trapped particles that are energetic enough to penetrate the walls and intervening 
material of the payload and be counted regardless of the counter axis orientation. 

The plots of counting rate versus X that are included herein are taken from 
only one of the passes analyzed. They are intended as a representative sample of the 
results obtained from this study. The plots for the other passes will be published 
separately. 

B. ERROR ANALYSIS 

Having obtained values of counting rate for various values of X, the accuracy 
of this determination must be investigated. Three sources of error are present in 
this analysis: error in reading the counting rate, error in establishing the time corre- 
sponding to a particular counting rate, and error in determining X as a function of 
time. 


The error in obtaining the counting rate becomes significant when the counting 
rate becomes high. Since the data are represented by a series of square waves whose 
periods are the time required to collect 2048 counts , a counting rate on the order of 
6000 counts sec”" 1 results in this square wave having a half period of approximately . 16 
second. The original data was recorded on strip charts with a paper speed of approxi- 
mately . 7 cm/sec. Hence, the subcarrier oscillator switches frequency every . 16 
second which corresponds to approximately 1 mm on the recorded strip chart. The 
location of these switches are probably not determined to better than + . 1 mm by the 
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FIGURE 29. COMPARISON OF OBSERVED COUNTING RATES WITH THE CALCULATED 
A. DURING THE INTERVAL 0500:00 to 0500:30 UT ON AUGUST 30. 










FIGURE 31. COMPARISON OF OBSERVED COUNTING RATES WITH THE CALCULATED 
X DURING THE INTERVAL 0510:00 TO 0510:30 UT ON AUGUST 30. THE 
PERIOD OF HIGH COUNTING RATE IS DUE TO A PASSAGE OF THE 
SATELLITE THROUGH AN ARGUS SHELL 
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FIGURE 35. COUNTING RATE VERSUS A FOR 0501: 30 TO 0502:00 UT, 
AUGUST 30, 1958 
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FIGURE 37. COUNTING RATE VERSUS A FOR 0510:07 TO 0510-25 UT 
AUGUST 30, 1958 



method used; hence, for high counting rates, the observed value may be in error by 
+ 10 per cent. Since this error is random and not cumulative, smoothing techniques 
may be used. This error could be reduced by playing back the original magnetic 
tapes and recording the data at a higher paper speed. A speed of 2 cm/sec should 
reduce this error to about + 3 per cent. 

The second source of error lies in establishing the times at which the counting 
rates are determined. Some variation in the paper speed at which the strip charts 
were, recorded was noted. To minimize this error , new reference times were taken 
every 15 seconds along the chart. It is believed that the variation in paper speed dur- 
ing a 1 5-second interval is not significant. However, these reference times cannot 
be located relative to Universal Time more accurately than to the nearest . 1 second, 
and at times when WWV is not clearly readable this precision cannot be guaranteed. 
During times when tumble modulations are present in the counting rate data, the 
counter axis may be moving relative to B at the rate of 60 deg/sec. A . 1-second error 
in time results in an error in X of about six degrees. 

The assumption that the observed counting rate corresponds to the midpoint 
between the two successive switches in the subcarrier oscillator frequency may also 
lead to a time error, particularly if the counting rate is changing rapidly. This error 
will be much less than the time interval between the two switches; therefore, for high 
counting rates this error will be less than . 1 second in time. Fow low counting rates, 
this error becomes greater, but low counting rates generally occur at small values 
of X. From the plots of counting rate versus X it was seen that counting rate is almost 
independent of X for 0° < \ < 50°; hence, this error may be disregarded. 

The final source of error lies in computing X as a function of time. The ac- 
curacy of this calculation obviously depends on the uncertainties in parameters <fi , ip, 
y, 6, bj, b 2 , and b 3 . The accuracy of <p has already been stressed. This is the most 
important parameter since cumulative errors may arise from this term. If c p is 
accurate to + . 0001 rad/sec, the resultant error in X will be less than four degrees 
after ten minutes. 

Since ip < < <p , the accuracy requirements on ip are not nearly as severe. In 
the cases where bj is small, all that is required of the ip term is to produce the long 
periods of nearly constant counting rates observed in the data. Since ip and 6 were 
obtained by fitting the observed counting rate data over the interval of interest with 
the condition that ip is constant, errors in these values will be neglected. It can be 
said that the requirement that ip be constant was consistent with the observed periods 
of counting rates where tumble modulations were not observed. 

The phase angle y involves the tumble axis L which has a fairly large uncer- 
tainty. However, the geometry of the L, U, and B during the_passes of interest is 
such that y is quite insensitive to errors in the orientation of L. The change in y re- 
sulting from a ten-degree variation in the right ascension and declination of L is on the 
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order of one degree, and will be neglected. An error in y results essentially in a time 
displacement of X. 

The uncertainty in L will also result in uncertainties in the bj, b 2 , and b 3 coefi- 
cients. Since t = 0 may be taken arbitrarily, b 3 = 0, b} = cos 0, and b 2 = sin /3 where /3 
is the angle B makes with L. In the passes examined, by chance, /3 was nearly 7r/2. 
Hence, when sin {ip t + 6) is nearly unity, X » /3, and the uncertainty in X is the same 
as in or about ten degrees. On the other hand, when cos {ip t + 6) is nearly unity, X 
is little affected by /3 when sin (<£ t + y) is small. When sin (<£ t + y) is nearly unity, 

X » 7r/2 - p. Hence, during the tumble modulations, X is well determined near -k/2. 

The minimum value that X reaches is subject to the same uncertainties involved in 
the L. Again, the fact that the radiation counting rate is insensitive to X for angles 
less than 50 degrees minimizes the effect of this error. 

C. CONCLUSIONS 

It has been demonstrated that the determination of the time dependence of 
the plane of tumble of Explorer IV has been successful in explaining the observed 
temperature variations, the observed drag accelerations, and the radiation counting 
rate modulations observed by the directional scintillation detector. Furthermore, 
counting rate as a function of angle that the scintillation counter axis makes with the . 
magnetic 'field has been determined at various points in the geomagnetic field. It has 
been shown that the counting rate is a maximum when the detector is perpendicular to 
the magnetic field and that the counting rate falls to a uniform value at angles less thar 
about 40 degrees through 50 degrees to the field. From the angular dependence of the 
counting rate, the angular distribution of the trapped particle flux may be obtained 
which will permit calculation of the intensity at each subsequent mirror point along 
the line of force. 

Two interesting effects were also observed during the course of this study. The 
first, which was originally observed by Fields, was that the roll required to explain 
the periodicity of the occurrences of intervals where the counting rate is not modulated 
by the tumbling motion appeared to vary froni day to day. The results of these obser- 
vations, shown in Figure 4, indicate that the roll actually increases by almost an order 
of magnitude on August 28. The roll rate increase observed on this date represents an 
angular momentum increase of 10 5 gm cm 2 sec -1 . It is not possible to determine 
whether total angular momentum is conserved and whether the increase in longitudinal 
angular momentum results from a decrease in transverse or tumble angular momen- 
tum, or if external torques are acting about the symmetry axis. A gain of 10 5 gm cm 2 
sec -1 about the longitudinal axis would result in a decrease in <p of approximately 
. 2 x 10 - 4 rad/sec if total angular momentum is to be conserved. This value is less 
than the uncertainty in <p. If external torques are responsible for the roll rate increase, 
a continuous torque on the order of one dyne-cm would be required. 
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The second interesting and somewhat surprising result is that the angular 
momentum vector does indeed appear to change orientation by as much as 10 degrees 
per day in an inertial reference frame. It has been recognized by Roberson [Ref. 52] , 
Doolin [Ref. 53] , and others that the gradient of the earth’s gravitational field can 
produce torques on a satellite which may result in a precession of the spin axis. An 
order of magnitude calculation of these gravity torques, using Roberson’s results, was 
made. The torques were time averaged over a tumble cycle and then time averaged 
over an orbital period. The resulting average torque was about 16 cyne-cm assuming 
an orientation of the tumble axis that would produce the maximum effect. Since the 
moment of inertia about the transverse axis is about 5 x 10 7 gm cm 2 and the tumble rate 
is approximately 1 rad/sec, the angular momentum is 5 x 10 7 gm cm 2 sec” 1 . A change 
in orientation of 10 degrees or . 17 rad represents an angular momentum change of 
. 85 x 10 7 gm cm 2 sec” 1 . For this change to occur in one day or 86, 400 seconds re- 
quires a continuous torque on the order of 10 2 dyne-cm. This exceeds the gravitational 
torques by a factor of six. 

Other effects that may give rise to precessing torques were examined cursorily. 
Aerodynamic forces near perigee will produce precessing torques due to the displace- 
ment of the center of lateral area from the center of gravity. The drag force at perigee, 
assuming an area of 3 x 10 3 cm 2 , an atmospheric density of 10" 13 gm cm -3 , a velocity 
of 8 x 10 5 cm sec^ 1 , and a drag coefficient of two, is about 200 dynes. Assuming a 
ten cm displacement between the center of area and center of gravity, the resultant 
torque is on the order of 2000 dyne-cm. However, this torque is only effective for the 
small percentage of time during an orbital revolution when the satellite is near perigee. 
Furthermore, after the satellite has begun tumbling in a plane, it would appear that 
this torque would average to zero over a tumble cycle. 

Interactions of a magnetic moment of the satellite with the geomagnetic field 
may give rise to precessing torques. If this magnetic moment is due to eddy currents 
induced by the tumbling motion of the satellite, the torques produced must be dissipative 
in nature. It has been established that the tumble rate changes by only 2 x 10” 8 rad 
see” 2 which corresponds to torques on the order of one dyne-cm. Therefore, the eddy 
currents appear to be much too small to produce the observed precessions. 

No measurements have been made on Explorer IV to determine whether or not 
it had a magnetic moment due to permanent magnetism, induced magnetism, or internal 
circuitry. Such measurements could be made on a spare payload and may or may not 
provide an explanation for the torques that appear to be acting on the satellite. 


63 



APPENDIX i 


The conservation of magnetic moment may be shown in the following manner. 
Starting with equation 10, 

d v 

m “dT” = B ’ 

multiplication by v J( gives 


_d 

dt 



dB d z 

d z d t 





hence, ju is constant with time. 
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APPENDIX 2 


Equation 13, 


<f> cos 0 


h 

h 


<^ 3 > 


may also be derived by direct geometrical argument. Let the total angular momentum 
L be along the z-axis. Let the body be rolling about I 3 with angular velocity o> 3 and 
precessing about L with the precession rate <p such that I 3 makes angle 0 with IT. 

Since angular velocities may be resolved into components and since I* is perpendicular 
to I 3 , the instantaneous angular momentum about Ij is I 2 <fi sin 0. The angular momenta 
about Ij and I 3 may be resolved into components along IT such that 


Ii <fi sin 


2 


0 + I 3 co 3 cos 0 


L, 


and components perpendicular to L such that 


Ij 0 sin 0 cos 0 - I 3 u> 3 sin 0 = 0. 
From the latte.r equation 

I3 ^3 

<p cos 0 = — ; • 
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"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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